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I. INTRODUCTION 

Rapidly increasing interest in the left-handed materi- 
als (LHM) started after Pendry et. at predicted that 
certain man-made composite structure could possess, in 
a given frequency interval, a negative effective magnetic 
permeability ficS^. Combination of such a structure with 
negative effective permittivity medium - for instance the 
regular array of thin metallic wires^"*" - enabled the con- 
struction of meta-materials with both effective permit- 
tivity and permittivity negative. This was confirmed by 
experiments^'^. 

Structures with negative permittivity and permittivity 
were named "left-handed" by Veselago^° over 30 years 
ago to emphasize the fact that the intensity of the electric 
field E, the magnetic intensity H and the wave vector k 
are related by a left-handed rule.^^ This can be easily seen 
by writing Maxwell's equation for a plane monochromatic 
wave: kx E = and k x H = -^E Once e and u, 

c c ^ 

are both positive, then E, H and k form a right set of 
vectors. In the case of negative e and /i, however, these 
three vectors form a left set of vectors. 

In his pioneering work, Veselago described the physical 
properties of LH systems: Firstly, the direction of the 
energy flow, which is given by the Poynting vector 



E X H 



(1) 



does not depend on the sign of the permittivity and 
permeability of the medium. Then, the vectors S and 
k are parallel (anti-parallel) in the right-handed (left- 
handed) medium, respectively. Consequently, the phase 
and group velocity of an electromagnetic wave propagate 
in opposite directions in the left-handed material. This 
gives rise to a number of novel physical phenomena, as 
were discussed already by Veselago. For instance, the 
Doppler effect and the Cherenkov effect are reversed in 
the LHMio. 

If both e and fj, are negative, then also the refraction 
index n is negative^*^'^^. This means the negative refrac- 
tion of the electro magnetic wave passing through the 
boundary of two materials, one with positive and the 
second with negative n (negative Snell's law). Obser- 
vation of negative Snell's law, reported experimentally^'^ 
and later in-'^^, is today a subject of rather controversially 
debate"'^^"^* . Analytical arguments of the sign of the re- 
fraction index were presented in^^. Numerically, nega- 
tive phase velocity was observed in FDTD simulations^^. 



Negative refraction index was calculated from the trans- 
mission and reflection data^°. Finally, negative refrac- 
tion on the wedge experiment was demonstrated also by 
FDTD simulations^!. 

Negative refraction allows the fabrication of flat lens^^. 
Maybe the most challenging property of the left-handed 
medium is its ability to enhance the evanescent modes^^ . 
Therefore flat lens, constructed from left-handed material 
with e = /i = — 1 could in principle work as perfect lens^^ 
in the sense that it can reconstruct an object without any 
diffraction error. 

The existence of the perfect lens seems to be in 
contradiction with fundamental physical laws, as was 
discussed in a series of papers^^'^^'^''. Nevertheless, 
more detailed physical considerations^^'^^'^^'^^'^^ not 
only showed that the construction of "almost perfect" 
lens is indeed possible, but brought some more insight 
into this phenomena^°~^^'^'''^^'^^'^^. 

As Veselago also discussed in his pioneering paper, the 
permittivity and the permeability of the left-handed ma- 
terial must depend on the frequency of the EM field, 
otherwise the energy density'^^ 
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would be negative for negative e' and /i' (real part of 
the permittivity and permeability). Then, according to 
Kramers- Kronig relations, the imaginary part of the per- 
mittivity (e") and of the permeability (^") are non-zero 
in the LH materials. Transmission losses are therefore 
unavoidable in any LH structure. Theoretical estimation 
of losses is rather difficult problem, and led even to the 
conclusion that LH materials are not transparent-'^'^. For- 
tunately, recent experiments^*'"^^ confirmed the more op- 
timistic theoretical expectation'^^, that the losses in the 
LH structures might be as small as in conventional RH 
materials. 

The number of papers about left-handed materials in- 
creased dramatically last year. The present paper is not 
the first review about left-handed materials. For recent 
reviews, see^" or the papers of Pendry"*!. Here we present 
typical structures of the left-handed materials (Sect. II), 
discuss a numerical method of simulation of the propaga- 
tion of EM waves based on the transfer matrix (Sect. Ill) , 
and present some recent results obtained by this method 
(Sects. IV, V). We discuss in Sect. V how the transmis- 
sion depends on various structural and material param- 
eters of LH structure. The method of calculation of the 
refractive index and of the effective permittivity and per- 
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meability is presented and applied to the LH structure. 
An unambiguous proof of the negative refraction index is 
given and the effective permittivity and permeabihty are 
calculated in Sects. VI, and VII. The obtained data for 
the permittivity are rather counter-intuitive and require 
some physical interpretation. Finally, in the Section VIII 
we discuss some new directions of the development of 
both theory and experiments. 
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FIG. 1. Left: Structure of the split ring resonator (SRR). 
The SRR consists of two splitted metallic "rings" . The SRR 
is characterized by the size w, width of the rings d, and two 
gaps: g and c. The external magnetic field induces an electric 
current in both rings^ . The shape of the SRR (square or circu- 
lar) is not crucial for the existence of the magnetic resonance. 
Right: Structure of the unit cell of the left-handed material. 
Each unit cell contains the split ring resonator located on 
the dielectric board, and one wire. Left-handed structure is 
created by regular lattices of unit cells. The EM wave prop- 
agates along the z direction. Periodic boundary conditions 
are considered in the x and the y direction, which assures the 
periodic distribution of the EM field. 



II. STRUCTURE 

LHM materials are by definition composites, whose 
properties are not determined by the fundamental physi- 
cal properties of their constituents but by the shape and 
distribution of specific patterns included in them. The 
route of the construction of the of LH structure consists 
from three steps: 

Firstly, the split ring resonators (SRR) (see fig. 1 for 
the structure of SRR) was predicted to exhibit the reso- 
nant magnetic response to the EM wave, polarized with 
H parallel to the axis of the SRR. Then, the periodic 
array of SRR is characterized^ by the effective magnetic 
permeability 



Mcfr(/) = 1 



(3) 



In (3), Vm is the resonance frequency which depends 
on the structure of the SRR (fig. 1) as (27r^'m)^ = 
3L2;Cj^igijj/[7rln(2c/d)r^]. F is the filling factor of the 
SRR within one unit cell and 7 is the damping factor 
27r7 — 2Lxp/r, where p is the resistivity of the metal. 

Formula (3) assures that the real part of is nega- 
tive at an interval around the resonance frequency. If 



an array of SRR is combined with a medium with nega- 
tive real part of the permittivity, the resulting structure 
would possess negative effective refraction index in the 
resonance frequency interval /S.v^^ . 

The best candidate for the negative permittivity 
medium is a regular lattice of thin metallic wires, which 
acts as a high pass filter for the EM wave polarized with 
E parallel to the wires. Such an array exhibits negative 
effective permittivity 



2,4,6 



eeff(i^) = 1 



with 



the 
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plasma 



frequency = Cjjgjj^/ (27ra^ ln(a/r))^. Sarychev and Sha- 
laev derived another expression for the plasma frequency, 
= Cught/(™^'C) with C = 21n(a/V2r)+7r/2 - 3^. 
Apart from tiny differences in both formulas, the two 
theories are equivalent^ and predict that effective per- 
mittivity is negative for v < Vp. 

By combining both the above structures, a left-handed 
structure can be created. This was done for the first 
time in the experiments of Smith et al.^ . Left-handed 
material is a periodic structure. A typical unit cell of 
the left-handed structure is shown in fig. 1. Each unit 
cell contains a metallic wire and one split ring resonator 
(SRR), deposited on the dielectric board. 
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FIG. 2. Transmission of the EM wave, polarized with E\\y 
and II through a periodic array of split ring resonators, 
wires, and of both SRR and wires. 



Fig. 2 shows the transmission of the EM waves through 
the left-handed structure discussed above. The transmis- 
sion through the array of the SRR is close to unity for 
all the frequencies outside the resonance interval (8.5-11 
GHz in this particular case) and decreases to -120 dB in 
this interval, because peS is negative (Eq. 3). The trans- 
mission of the array of metallic wires is very small for all 
frequencies below the plasma frequency (which is ^ 20 
GHz in this case), because Coff is negative (Eq. 4). The 
structure created by the combination of an array of SRR 
and wires exhibits high transmission T ^ 1 within the 
resonance interval, where both e^s and p^s are negative. 
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For frequencies outside the resonance interval, the prod- 
uct /Ueffeoff is negative. The transmission decays with the 
system length, and is only ~ — 120 dB in the example of 
fig. 2. Experimental analysis of the transmission of all 
the three structures was performed by Smith et al.^ . 

We want to obtain a resonance frequency Vm ~ 10 
GHz. This requires the size of the unit cell to be 3-5 
mm. The wavelength of the EM wave with frequency 
~ 10 GHz is w 4 cm, and exceeds by a factor of 10 the 
structural details of the left-handed materials. We can 
therefore consider the left-handed material as macroscop- 
ically homogeneous. This is the main difference between 
the left-handed structures and the "classical" photonic 
band gap (PEG) materials, in which the wave length is 
comparable with the lattice period. 

It is important to note that the structure described in 
fig. 1 is strongly anisotropic. For frequencies inside the 
resonance interval, the effective Ceff and are negative 
only for EM field with H \\ x and E || y. The left-handed 
properties appear only when a properly polarized EM 
wave propagates in the z direction. The structure in fig. 
1 is therefore effectively one-dimensional. Any EM waves, 
attempting to propagate either along the x or along the 
y direction would decay exponentially since the corre- 
sponding product eeff/ieff is negative. This structure is 
therefore not suitable for the realization of the perfect 
lens. To test the negative Snell's law experimentally, a 
wedge type of experiment must be considered^^ (fig. 3) 
in which the angle of refraction is measured outside^ the 
left-handed medium (in air)^'^'^'*. 



Two dimensional structures have also been con- 
structed. For instance, the anisotropy in the x — z plane 
is removed if each unit cell contains two SRR located 
in two perpendicular planes^'"'^^. No three dimensional 
structures have been experimentally prepared yet. 




FIG. 4. Two positions of the metallic wire in the unit cell. 
Left: wire is deposited on the opposite side of the dielec- 
tric board®'^"*. Right: wire is located along the split-ring 
resonator*^ . 



It is also worth mentioning that some other one dimen- 
sional structures were prepared, in which the wires were 
deposited on the same side of the dielectric board with 
the SRRs. The wires could be located either on the oppo- 
site side of the dielectric board, as it was done in9'i3,38^ 
or put next to the SRRs on the same side of board**^ 
as shown in the right panel of fig. 4. More complicated 
one-dimensional structure was suggested by Ziolkowski*'^. 
Recently, Marques et al. found left-handed behavior in 
an array of SRR located inside a metallic wave guide*'*. 




FIG. 3. Refraction experiment on the Left-handed 
material The incident EM wave propagates from the left 
and hits perpendicularly the left boundary of the wedge. The 
angle of refraction is measured when the EM wave passes the 
right boundary of the inspected material and propagates for 
some time in the air. Two possible directions of the propaga- 
tion of the refracted wave are shown: the right refraction for 
the conventional right-handed (RH) material, and left refrac- 
tion for the left-handed material. This experimental design 
enables to use also strongly anisotropic one-dimensional LH 
samples, since the angle of refraction is measured outside the 
sample. 



III. NUMERICAL SIMULATION 

Various numerical algorithms were used to simulate 
the propagation of EM waves through the LH structure. 
We concentrate on the transfer matrix algorithm, devel- 
oped in a series of papers by Pendry and co-workers*^. 
The transfer matrix algorithm enables us to calculate 
the transmission, refiection, and absorption as a func- 
tion of frcquency*^'*^. Others use commercial software: 
either Microwave studio''*^''^* or MAFIA^'*°, to estimate 
the position of the resonance frequency interval. Time- 
dependent analysis, using various forms of FDTD algo- 
rithms are also used**'*^'*^'*^'^". 

In the transfer matrix algorithm, we attach in the z 
direction, along which EM wave propagates, two semi- 
infinite ideal leads with e = 1 and /i = 1. The length 
of the system varies from 1 to 300 unit cells. Periodic 
boundary conditions along the x and y directions are 
used. This makes the system effectively infinite in the 
transverse directions, and enables us to restrict the sim- 
ulated structure to only one unit cell in the transverse 
directions. 

A typical size of the unit cell is 3.66 mm. Because 
of numerical problems, we are not able to treat very 
thin metallic structures. While in experiments the thick- 
ness of the SRRs is usually 17 ^m, the thickness used in 
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the numerical calculations is determined by the minimal 
mesh discretization, which is usually « 0.33 mm. In spite 
of this constrain, the numerical data are in qualitative 
agreement with the experimental results. This indicates 
that the thickness of the SRR is not a crutial parameter, 
unless it decreases below or is comparable with the skin 
depth 5. As (5 ~ 0.7/im at GHz frequencies of interest, 
we are far from this limitation. We will discuss the role 
of discretization in Section IV. 
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IV. TRANSMISSION 

As discussed in Section II the polarization of the EM 
waves is crucial for the observation of the LH properties. 
The electric field E must be parallel to the wires, and 
the magnetic field H must be parallel to the axis of the 
SRRs. In the numerical simulations, we treat simultane- 
ously both polarizations, E\\x as well as E \\ y. Due to 
the non-homogeneity of the structure, these polarizations 
are not separated: there is always non-zero transmission 
txy from the x to y polarized wave. As we will see later, 
this effect is responsible for some unexpected phenomena. 
At present, we keep in mind that they must be included 
into the formula for absorption 



1 l^xa;! Ka;?/! ^ 



xy\ 



(5) 



and in the equivalent relation for Ay. 

Figure 5 shows typical data for the transmission in the 
resonance frequency region. A resonance frequency inter- 
val, in which the transmission increases by many orders 
of magnitude is clearly visible. Of course, high transmis- 
sion does not guarantee negative refraction index. The 
sign of n must be obtained by other methods, which will 
be described in Sect. VI. 

In contrast to the original experimental data, numer- 
ical data show very high transmission, indicating that 
LH structures could be as transparent as the "classical" 
right-handed ones. This is surprising, because due to the 
dispersion, high losses are expected. Fig. 6 shows the 
transmission as the function of the system length, for 
three different frequencies inside the resonance frequency 
interval. 
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FIG. 5. Transmission power T through the Left-handed 
meta material of various lengths. Left: frequency dependence, 
right: length dependence. 
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FIG. 6. System length dependence of the real part of the 
transmission through the left-handed structure for three var- 
ious frequencies within the left-handed band. The system 
length is given as the number of unit cells in the propagation 
direction. Note the different scale on the y axis for the three 
cases. 

Figure 5 shows also that the transmission never 
decreases below a certain limit. Due to the non- 
homogeneity of the structure there is a non-zero prob- 
ability tyx that the EM wave, polarized with E \\ y, is 
converted into the polarization E \\ x. The total trans- 
mission tyy consists therefore not only from the "unper- 



"yy 

turbed" contribution t'fy., but also from additional terms, 
which describe the conversion of the y-polarized wave 
into x-polarized and back to y-polarized: 

tyy{Q, L) = t,(°)(0, L) + V txy{Q, z%x{z, z')tyx{z'L) + ... 



(6) 



t^yy (0, L) decreases exponentially with the system length 
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L, while the second term, which represents the conversion 
of the y-polarized wave into a:-polarized wave and back, 
remains system-length independent, because txx{z,z') ~ 
1 for any distance \z — z'\. 
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FIG. 7. The transmission peak for various sizes of the unit 
cell. In contrast to the structure shown in Fig. 1, the metallic 
wire is located on the opposite side of the dielectric board. 
This enables to compress the width of the unit cell as shown 
in the far left panel. 



V. STRUCTURE AND TRANSMISSION 

Numerical simulations confirm the existence of the res- 
onance left-handed frequency interval. Before we proceed 
in the calculation of the effective system parameters, let 
us briefly discuss how the structure of the unit cell influ- 
ences the position and the width of the resonance inter- 
val. 
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FIG. 8. Top panel: Effective permittivity eoff for an ar- 
ray of metallic wires with radius 0.1 mm. The lattice con- 
stant (the distance between the wires) is 5 mm. The electric 
field is parallel to the wires. Note that the imaginary part 
of the effective permittivity is multiplied by a factor of 10. 
We present results of the analysis of the transfer matrix data 
for three different discretization of the unit cell. The lower 
panel shows the ratio \e" /e\ for an array of thin metallic wires 
100 X 100 ^m for two values of the metallic permittivity. The 
numerical results prove that electromagnetic losses decrease 
when the metallic conductance increases. 

The resonance frequency depends on the structural pa- 
rameters of the SRR and on the parameters of the unit 
cell^ . Qualitative agreement between the theoretical for- 
mulas and the numerical results was obtained'*^'^^. An- 
other problem is the dependence of the resonance interval 
on the size and shape of the unit cell. As an example, 
we show in Fig. 7 how the transmission through the var- 
ious LH structures depends on the width/high/length of 
the unit cell. Notice that the width of the LH frequency 
interval increases substantially as shown in the far left 
panel of fig. 7. 

Another important parameter is the permittivity em 
of the metallic components. Within the first approxima- 
tion, we can consider both the SRR and the wires made 
from a perfect metal. Then both the conductivity a and 
the imaginary part of the metallic permittivity are 
infinite''^ . This option is often used in the simulations 
of the commercial software^*. In the transfer matrix al- 
gorithm, however, is finite, of the order of 10^. For 
copper, which is currently used in the experimentally fab- 
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ricate LH structures, ej^, « lO''^^. Some test calculations 
with higher values of e'^ gave us almost the same result, 
indicating that ej^ ~ 10^ is already sufficient to simulate 
realistic materials'^^. 

As was mentioned in Sect. I, the left-handed struc- 
ture must be dispersive. Dispersion requires a non-zero 
imaginary part for the permittivity and permeability. 
Therefore, transmission losses can not be avoided in LH 
systems. This seems to be in agreement with the first 
experimental data: the transmission measured in the 
experiments^'^, was only of the order 10^'^. Although re- 
cent experiments reported the transmission very close to 
unity, there are still serious doubts in the literature^^'^^ 
about the possibility to create highly transmitted LH 
structures. Results of numerical simulations, as shown 
in fig. 5, however give that the transmission through the 
LH structure could be very high, of order of unity. 

The main argument against the expectation of high 
transmission in LH structures'^ is that the effective per- 
mittivity of the array of thin metallic wires is mostly 
imaginary than real and negative^^. This is a serious 
objection, because the thickness of the wires in fig. 1 
is 1 mm, which exceeds more than 10 times the real- 
istic parameters of the experimentally analyzed struc- 
tures. There is no chance to simulate very tiny metallic 
structures in the transfer matrix algorithm. Fortunately, 
Pendry et al. (formula (36) of^) had shown that the prob- 
lem of the thickness of the wires could be avoided by 
simultaneously re-scaling the metallic permittivity. An 
increase of the wire radius by factor a could be compen- 
sated by a decrease of the permittivity of the metal by 
factor of . As the metallic permittivity used in^° is 
« 100 times lower than the realistic permittivity of Cu, 
we assume that the size of the wires (1 mm) corresponds 
to Cu wires of thickness 0.1 mm, which is close to size 
used in experiments (0.2 x 0.017 mm). 




FIG. 9. The dependence of the transmission peak on the 
imaginary part of the dielectric permittivity of the dielectric 
board. Note that standard dielectrics have e" « 10"'^. 

In fig. 8 we present results for the effective permittivity 
of an array of thin metallic wires. The left panel confirms 
that the transfer matrix algorithm provides us with real- 



istic data for the transmission and reflection. Three dif- 
ferent discretizations of the unit cells were used. For each 
of them the effective permittivity was calculated by the 
method explained in Sects. VI, VII. The right panel of 
figure 8 shows that transmission losses are small in real- 
istic LH structures'^. Losses in the metallic components 
of the left-handed structure are therefore not responsible 
for the low transmission, measured in the experiments*'^. 
The role of the material properties of the dielectric board 
on which SRR is deposited was also studied'^^. As is 
shown in Fig. 9, very small imaginary part of eeoard 
causes a rapid decrease of the transmission. This surpris- 
ing result seems to be in agreement with experiments^*. 
Much higher transmission was obtained in cases with ex- 
tremely small losses in the dielectric board^*. 
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FIG. 10. To panel shows n'kL as a function of the system 
length L for the EM wave polarized with E parallel to the 
wires. The slope is negative which confirms that n' is neg- 
ative. Bottom: n'kL as a function of the system length L 
for the EM wave polarized E \\ x. The slope is positive which 
confirms that Re n' is positive. A linear fit gives that n' = 1.1. 
There is almost no interaction of the LH structure with the 
EM wave. 



VI. EFFECTIVE INDEX OF REFRACTION 

As was discussed above, the structural inhomogeneities 
of the LH materials are approximately ten times smaller 
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than the wavelength of the EM wave. It is therefore pos- 
sible, within a first approximation, to consider the slab of 
the LH material as an homogeneous material. Then we 
can and use the textbook formulas for the transmission 
te^'^^ and reflection r for the homogeneous slab: 



and 



cos(nfci) ~ ^ f + ) sin(n/;;Z/) 



(7) 



(8) 



Here, z and n are the effective impedance and the re- 
fraction index, respectively, k is the wave vector of the 
incident EM wave in vacuum, and L is the length of the 
LH slab. We only consider perpendicular incident waves, 
so that only the z components of the effective parameters 
are important. To simplify the calculations, we neglect 
also the conversion of the polarized EM wave into an- 
other polarization, discussed in the Section IV. More 
accurate analysis should treat both t and r as 2 x 2 ma- 
trices. Here we assume that the off-diagonal elements of 
these matrices are negligible: 



(9) 



In the present analysis, we use the numerical data for 
the transmission and the reflection obtained by the trans- 
fer matrix simulation. The expressions for the transmis- 
sion and the reflection can be inverted as 



2 = ±4 



(1 - r)2 - t2 



cos(nfcL) = X = — (1 - + 1^) 



(10) 



(11) 



The sign of z is determined by the condition 

z' >0 (12) 

which determines z unambiguously. The obtained data 
for z enables us also to check the assumption of the ho- 
mogeneity of the system. We indeed found that z is in- 
dependent of the length of the system L. 

The second relation, (11), is more difficult to invert 
since cos~^ is not an unambiguous hmction. One set 
of physically acceptable solutions is determined by the 
requirement 



n" > 



(13) 



which assures that the material is passive. The real part 
of the refraction index, n', however, suffers from the un- 
ambiguity of 27rm/(fci) (m is an integer). To avoid this 
ambiguity, data for various system length L were used. 
As n characterizes the material property of the system. 



it is L independent. Using also the requirement that 
n should be a continuous function of the frequency, the 
proper solution of (11) was found and the resonance fre- 
quency interval, in which n' is negative was identified^". 
Here, we use another method for the calculation of n and 
z: Equation (11) can be written as 



6"""'=^ [cos(n'fcL) + i sin(n'fcL)] = Y = X± y/l-X'^. 

(14) 

Relation (14) enables us to find unambiguously both the 
real and the imaginary part of the refraction index from 
the linear fit of n"kL and n'kL vs the system length L. 
The requirement (13) determines the sign in the r.h.s. of 
Eq. 14, because \Y\ < 1. Then, the linear fit of n'kL vs 
L determines unambiguously the real part of n'. 

Figure 10 shows the L dependence of n'kL for three fre- 
quencies inside the resonance interval (transmission data 
for these frequencies are presented in fig. 6). The numer- 
ical data proves that the real part of the refraction index, 
n' . is indeed negative in the resonance interval. For com- 
parison, we present also n'kL vs L for the x polarized 
wave outside the resonance interval. As expected, the 
slope is positive and gives that n' = 1.13, which is close 
to unity. 

Besides the sign of the real part of the refraction in- 
dex, the value of the imaginary part of n, n" is important, 
since it determines the absorption of the EM waves in- 
side the sample. Fortunately, n" is very small, it is only 
of the order of 10~^ inside the resonance interval. As 
is shown in fig. 6, quite good transmission was numer- 
ically obtained also for samples with length of 300 unit 
cells (which corresponds to a length of the system 1.1 
m). This result is very encouraging and indicates that 
left-handed structures could be as transparent as right- 
handed materials. 

There are two constrains in the present method, (i) 
While the above method works very well in the right 
side of the resonance interval, we had problems to es- 
timate n in the neighborhood of the left border of the 
resonance region, where n' is very large and negative. 
This is, however, not surprising since in this frequency 
region the wavelength of the propagating wave becomes 
comparable with the size of the unit cell, so that the effec- 
tive parameters have no physical meaning, (ii) Outside 
the resonance interval, we have serious problems to re- 
cover proper values of n. This is due to the conversion of 
the X polarized wave into a y polarized. As a result, we 
do not have enough numerical data for obtaining the L- 
dependence of the transmission t. As it is shown in Fig. 
5, only data for 2-3 unit cells are representative for the 
E\\y wave. The condition (9) is not any more fulfilled, 
and the present theory must be generalized as discussed 
above. 
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FIG. 11. Refraction index n, impedance z, permittivity e 
and permeability /i as a function of frequency in the resonance 
interval. Close to the resonance frequency v = 9.8GHz, there 
are problems to estimate n and z, because the wave length of 
the propagated EM wave becomes comparable with the size of 
the unit cell of the left-handed structure. Dashed area shows 
the resonance interval. 



VII. EFFECTIVE PERMITTIVITY AND 
PERMEABILITY 



Once n and z are obtained 



Ccff = n/z 



and 



(15) 



one can estimate unambiguously the effective permittiv- 
ity and permeability of the investigated structure. Re- 
sults are shown in fig. 11. Frequency dependence of 
the effective refraction index and e'^g, /Lt^g is also pre- 
sented in^*' . Our numerical results confirm the resonance 
behavior of the fi'^g, in agreement with theoretical pre- 
dictions of Pendry et . In addition to these expected 
results, we also found rather strong electric response of 
the SRR, which manifests itself as the "anti-resonant" 
behavior of e^^ in the resonance frequency interval and 
by the decrease of the plasma frequency (fig. 4 of ref.^°, 
see also^^). 

The real parts of the permittivity and permeability, 
e^g and /i^jj are negative in the resonance interval, as ex- 
pected. Surprisingly, the imaginary part of the effective 
permittivity, e"g was also found negative. This seems to 
contradict our physical intuition, since, foUowing^^, Sect. 
80, the electro-magnetic losses are given by 



Q=l.Jdu;co [e"\Ef + ^,"\Hf] 



(16) 



In passive materials, we require Q > 0. This is trivially 
satisfied if we require that both e" and /i" are positive. 
Negative e"ff was obtained also in other structures'^. 

More detailed analysis is also necessary for the energy 
of the EM wave inside the left-handed material. The for- 
mula for the total energy, given by Eq. 2, does not assure 
the positiveness of the energy because the first term in 
Eq. (2) is negative in the left part of the resonance in- 
terval. Indeed, as shown in fig. 11, both e' and de' /dcu 
are negative, so that 



djuje') 
duj 



< 0. 



(17) 



Although formulas (2,16) are not the most general for- 
mulas for the energy of the EM field'^ , we show that both 
e'^fi < and relation (17) are consistent with (16) and (2), 
respectively. We use the definition of the impedance^^. 
Sect. 83, 



(18) 



8 



Then, with the help of relations (15), we re- write the 
relation (16) into the form 

Q = J- / diouj\H\'^ X 2n"{Lj)z'{Lj) (19) 

ZTT J 

which assures that Q is positive thanks to the conditions 

that z' > and n" > (12,13). There is therefore no 
physical constrain to the sign of the imaginary part of 
the permittivity and permeability. 

With the help of the relation (18) Eq. 2 can be rewrit- 
ten into the form 



U 



27r/ 



(20) 



Using the numerical data for the impedance and for real 
part of the permittivity, we checked that indeed the ex- 
pression in the bracket of the r.h.s. of (20) is always 
positive. 



VIII. FURTHER DEVELOPMENT, UNSOLVED 
PROBLEMS AND OPEN QUESTIONS 

We reviewed some recent experiments and numerical 
simulations on the transmission of the electromagnetic 
waves through left handed structures. For completeness, 
we note that recently, Notomi^^ has studied the light 
propagation in strongly modulated two dimensional pho- 
tonic crystals (PC). In these PC structures the permit- 
tivity is periodically modulated in space and is positive. 
The permeability is equal to unity. Such PC behaves as a 
material having an effective refractive index controllable 
by the band structure. For a certain frequency range it 
was found by FDTD simulations^^"^'' that rieg is nega- 
tive. It is important to examine if left-handed behavior 
can be observed in photonic crystals at optical frequen- 
cies. 

Negative refraction on the interface of a three dimen- 
sional PC structure has been observed experimentally by 
Kosaka et al.^^ and a negative refractive index associ- 
ated to that was reported. Large beam steaming has 
been observed in^^, that authors called "the supcrprism 
phenomena" . Similar unusual light propagation has been 
observed in one-dimensional and two dimensional refrac- 
tion gratings. Finally, a theoretical work^° has predicted 
a negative refraction index in photonic crystals. 

Studies of the left handed structures open a series of 
new challenging problems for theoreticians as well as for 
experimentalists. The complete understanding of the 
properties of left-handed structures requires the reevalu- 
ation of some "well known" facts of the electromagnetic 
theory. There is no formulas with negative permeabil- 
ity in classical textbooks of electromagnetism^^'^^. Ap- 
plication of the existing formulas to the analysis of left 
handed structures may lead to some strange results. The 
theory of EM field has to be reexamined assuming neg- 
ative /U and e. We need to understand completely the 



relationship between the real and the imaginary parts of 
the permittivity and the permeability. Kramers-Kronig 
relations should be valid, but nobody have verified them 
yet in the case of the left-handed structures. The main 
problem is that we need eeS and fieS in the entire range 
of frequencies, which is difficult to obtain numerically. 
Then, due to the anisotropy of the structure as well as 
the nonzero transmission txy in Eq. (7), Kramers-Kronig 
relations should be generalized. We do not believe that 
today's numerical data enables their verifications with 
sufficient accuracy. 

Both in the photonic crystals and LHM literature there 
is a lot of confusion about what is the correct definitions 
of the phase and group refractive index and what is their 
relations to negative refraction. In additions, it is in- 
structive to see how the LH behavior is related with the 
sign of the phase and group refractive indices for the PC 
system. The conditions of obtaining LH behavior in PC 
were recently examined in^^. It was demonstrated that 
the existence of negative refraction is neither a prerequi- 
site nor guarantees a negative effective refraction index 
and so LH behavior. Contrary, LH behavior can be seen 
only if phase refractive index riphasc is negative. Once 
^phase is negative, the product S.k is also negative, and 
the vectors k, E and H form a left handed set, as dis- 
cussed in the Introduction. 

Problems of causality arises also in connection with 
"well 

known" and accepted relations like d{n' {uj)uj)/duj > 0. 
It is evident, that this relation can not be valid in the 
vicinity of the left border of left-handed frequency inter- 
val, since both n' and dn! jdio are negative. The same 
problem arises for the real part of Ceff. Also relation (17) 
requires more exact and complete treatment. We need 
more general relations for the energy of the EM field^^ 
which incorporates all the allowed signs of the real and 
imaginary parts of the permittivity and permeability. 

Following problems that are currently discussed in lit- 
erature. The negative Snell's law requires the under- 
standing in more detail of the relationship between the 
Poynting vector, the group velocity, and the phase ve- 
locity. We believe, that there is no controversial in this 
phenomena Anisotropy of real left-handed materials 
inspires further development of super- focusing''^. 

We believe that further analysis, of what happens when 
EM wave crosses the boundary of the left-handed and 
right-handed systems, will bring more understanding of 
the negative refraction as well as perfect lensing. Nu- 
merical FDTD simulations of the transmission of the EM 
wave through the interface of the positive and negative 
refraction index'"'^ showed that the wave is trapped tem- 
porarily at the interface and after a long time the wave 
front moves eventually in the direction of negative refrac- 
tion. Computer simulations of the transmission through 
LH wedge^"'^ also confirm that EM wave spends some time 
on the boundary before the formation of the left-handed 
wave front. Formation of surface waves^"'"^^'^^ can ex- 
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plain "perfect lensing" without violating causality. Re- 
cent development, both numerical'''*' and experimental"''* 
indicate that perfect (although not absolutely perfect) 
lensing might be possible. We need also to understand 
some peculiar properties of the left handed structures due 
to its anisotropy''^ and bi-anisotropy^^. 

Last but not least, let us mention an attempt to 
find new left-handed structures, both in the traditional 
metallic left-handed structures'^ and in the photonic 
crystals-'^'^'*^''. 
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